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The article reviews recent results for the low energy physics of 
fast tunneling centers in metallic environments. For strong enough 
couplings to the environment these tunneling centers display an orbital 
Kondo effect and give rise to a non-Fermi-liquid behavior. This latter 
property is explained by establishing a mapping of the tunneling center 
model to the multichannel Kondo model via the renormalization group 
transformation combined with a 1/Nf expansion. The case of M- 
state systems, the role of the splittings and the present experimental 
situation are also discussed. 

^ i PACS numbers: 72.10.Fk, 72.15.Cz, 71.55.-i 
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Non-Fcrmi-liquid systems like some heavy fermion mater onefdimensional 
interacting electronsa or multichannel quantum impurity problemsuta have been 



I 

q . the subject of growing interest during the past few years. These systems have the 
(j common feature that they display nonanalytical power law singularities in various 
, physical-properties and that the usual quasiparticle picture of Landau's Fermi- liquid 
I ' theorytffl does not apply for them. 

From the systems mentioned above quantum impurity problems are of special 
interest. This is due to the fact that they provide simple toy models for the study 
Q^ ' of strongly correlated lattice models agd furthermore the techniques developed fqx 



them can be directly applied to finiteoa and infinite dimensional lattice models E 
These latters seem to be very close to the realistic-two or three-dimensional models 
, designed to describe heavy fermion compounds S3 Furthermore it has been sug- 
' gested that even the properties of hlsh-T c materials might be explained by means 
— ^ ' of multichannel Kondo-like models .Eil 

One of the most interesting quantum impurity problems is provided by tunneling 
centers (TC). These are formed by some heavy particle (an atom, a point defect 
I or a collective coordinate of a dislocation) moving in an effective potential with at 
*^ ' least two local minima. If the barrier between these minima is sufficiently large and 
Ch , the minima are close enough to each other then at low temperature the hopping via 
thermal activation is negligible and the atom can move only by tunneling from one 
minimum to another. At low enough temperatures the heavy particle stays in the 
lowest lying quantum states associated to the minima of the effective potential and 
the presence of the higher excited states is only reflected in the renormalization of 
the different coupling const ants characterizing the interaction of the TC with its 



03 surrounding (See Sec. VII E| ) 



The physics of an isolated TC is rather boring, but becomes interesting if the 
TLS is coupled to its environment. In this paper we restrict our considerations to 
TC's put into a metal. A realistic TC in a metal is coupled both to the acoustic 
phonons and the conduction electrons providing the low energy excitations of the 
environment. However, as also stressed by Prokof'ev,c3 since the phonon-TC cou- 
pling is proportional to the momentum transfer q and the density of the low-energy 
phonon modes scales as ~ ui 2 this interaction can be neglected with respect to the 
coupling to the low energy electron-hole excitations, which have a linear density of 
states and an energy independent coupling to the TC. |— ,|— , 
The TC is coupled to the conduction electrons by two different processes £3ii-3 
For a fixed position of the heavy particle the conduction electrons tend to build 
up a screening cloud around the tunneling particle. Since by a tunneling process 
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the heavy particle has to carry with itself this screening cloud consisting of an 
infinite number of electron-hole excitations, this interaction tends to localize the 
heavy particle reducing its tunneling rate at low temperatures. For slow tunneling 
centers displaying individual jumps this is the only important process-and it has 
been studied in great detail both experimentally and theoreticallyOESO If only 
this screening interaction is taken into account then the different vertices occurring 
in a perturbative expansion commute with each other. Therefore this model is 
sometimes also referred to as the commutative model. 

On the other hand, if the tunneling rate is large enough, then in addition to 
the sc reenin g electron assisted tunneling becomes also possible (noncommutative 
model) JHrEa In this process an electron is scattered on the heavy particle while 
it jumps from one minimum to another. The latter is essential!}; due to barrier 
fluctuations caused by the local electronic density fluctuationsE§Ej but, it can also 
be generated by the virtual hoppings to the excited states of the TC.c3 While the 
amplitude of this process is rather small, it results in the generation of a Kondo 
effect and drives the system away from the marginally stable commutative model 
at low energy scales. As the temperature decreases, the conduction ele ct l j oi ns and 
the heavy particle form a strongly correlated Kondo-like bound state. EJ'Hill For a 
two-level system (TLS), the important case of a TC with two potential minima, 
the energy of this ground state is approximately given by the Kondo energy 



where D is some bandwidth cutoff of the order of the Fermi energy, and v z and v x 
denote the dimensionless couplings characteristic to screening and assisted tunnel- 
ing, respectively. For an M-level system no such simple formula can be given, but 
the estimated maximal Kondo temperatures are of the same order of magnitude, 



In the considerations of the previous paragraph we have neglected the role of 
the splitting A of the nearly degenerate levels of the TC. While this splitting is 
usually small and it is also renormalized downward during the scaling procedure 
it generates a further low energy scale called the freezing temperature T*. Above 
this scale the levels of the TC can be considered as degenerate. On the other hand, 
below T* the splitting becomes relevant and the mlernal degrees of freedom leading 
to the Kondo effect are in most cases frozen out £3 

Usually one assumes that the interaction of heavy particle and the conduction 
electrons is independent of the spin of the electrons, thus the spin up and spin 
down conduction electron channels are scattered independently and exactly with 
the same amplitude by the heavy particle. While the spin is only present as a silent 
quantum number (flavor) it plays a crucial role. As discussed later, in the appro- 
priate temperature range T* < T < Tk the low energy properties of an M-state 
TC caft-he, described by the two channel SU(M) x SU(N f = 2) Coqblin-Schrieffer 
modelJl3'El'E3c3 a *** generalization of the multichannel Kondo model. Here the 
SU(M) symmetry is connected to the M positions of the TC, while the SU(2) 
symmetry is due to the spins of the conduction electrons. This model, in contrast 
to the usual SU (M) single channel model exhibiting a Fermi-liquid behavior ,0oE3 
belongs to the family of overcompensated spin models and displays non-Fermi- 
liquid properties due to the iVy-fold degeneracy in the 'flayer'., and gives rise to 
nonanalytical behavior of the different physical quantities.E2rEJt!3H In the special 
case of a TLS (M = 2) the equivalent model is just the two-channel Kondo model. 
While the usual SU(2) x SU{2) two-chamm| -Kga do model has been investigated by 
various techniques such as Bethe AnsatzJi§EMiconformal field theory,r 4 J large N 
expansions,!! numerical renormalization group,Ej path integral approach ,EjO and 
bosonazitaions,c3 the number of results obtained for the SU(M) x SU(N) model 
with M > 3 is still more restrictive. 

In the present review we mainly concentrate on the mapping of the noncommuta- 
tive TC model to the SU (M) x SU (N) model in the above described temperature 
range, but we also give a concise review of the recent experimental and theoretical 
results in the field. 

The paper is organized as follows: In Section O we shortly discuss the model 




(1.1) 
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introduced to describe a tunneling heavy particle. In Sec 
gin of the logarithmic corrections appearing in the p 
Anderson's 'poor man's scaling' to the TLS problem 



III we discuss the ori- 



urbation theory and apply 
Then in Sec. |V] we give a 
concise introduction to the multiplicative renormalizat ion group and we show how 
to obtain the estimation of the Kondo temperature (1.1). Section ^ is devoted 
to the development of a large Nf expansion in the 'flavor degeneracy', Nf for the 
TLS. We show that up to 1/Nj order the generalized 7V/-fiavor TLS model can be 
mapped to the _/V/-channel Kondo model. The generalization of the above map- 
ping to the case of M-state systems will be discussed in Sec. VL Some additional 
results not closely connected to the subject of this paper such as the path integral 
formalism applied to two-level systems, the influence of spin-flip sca tteri ng on the 
TLS, and the role of excited states will be shortly discussed in Sec. VII. We also 
discuss the present experimental situation and the poss ibility of the observation 
of the non- Fermi-liquid behavior, predicted in Sec. VIII. Finally, our concluding 
remarks are given in Section IX. 



II. THE TWO-STATE MODEL 



To investigate the behavior of two-level system (TLS) coupled to the conduction 
electrons we first have to construct a model that contains all the relevant ingredients 
of the realistic situation. In the present Section we restrict our considerations to 
the case of TLSs but they can easily be extended to the general case of M-level 
systems (MLSs). 

In the nature a variety of physical realizations of TLS's exist. However, to be 
specific, constructing the model we consider the case of amorphous metals being 
one of the most frequently studied systems containing TLSsE3 The construction of 
the model for other systems follows similar lines. In Fig[|.a we represent a situation 
when one or several neighboring atoms (ions) in the amorphous structure have two 
stabile positions close to each-other. This situation can be most easily described by 
some effective potentia -shown in Fig. ^.b, where R denotes the relevant variable in 
the configuration spaceE_3 and the two minima are associated to the two positions of 
the TLS in Fig. |l].a. Then one can associate two quantum states <f>i(R) and 4> r {R) 
of the atoms to the two minima of the potential well, which are localized in the 
left- and right-hand side, respectively. We assume, that the shape of the double 
well potential is such that there are two low-energy states with energies E\ ~ E2 '• 
two linear combinations of <fii(R) and <fr r (R), and the next one has a much higher 
energy: £3 — E\ ~ WDcbye- With this assumption the motion of the TLS at low 
temperatures will be constrained to the two lowest lying states and therefore the 
states having higher energies can be ignored. In an amorphous metal the energy 
distance A between the two lowest states can differ from TLS to TLS. 





Fig. a 



Fig. b 



FIG. 1. The TLS in an amorphous structure and the effective potential describing its 
motion. 



In the following considerations we restrict ourselves to a single TLS and an elec- 
tron band of the simplest shape. To obtain macroscopic quantities we shall treat 
each TLS as being independent of the other TLS's and in the end of the calculation 
we perform an averaging over the distribution of the individual TLS parameters. 
The separation of this last summation is reliable at low TLS concentration, i.e., 
when they are far from one another and do not interact. 

At low enough temperatures the Hamiltonian for a TLS could be described simply 
by a 2 x 2 matrix acting in the space of the states <fii(R) and <fi r (R). Nevertheless, 
for later convenience we rather describe the TLS motion in ±erms of an imaginary 
fermion (pseudofermion) , initially introduced by AbrikosovEj to make the second 
quantized formalism applicable for the Kondo problem. In this language the physi- 
cal states 4>i{R) and (f> r (R) of the TLS correspond to the one particle pseudofermion 
states b^\0) and 10), respectively, b^ and 6+ being the pseudofermion creation 
operators. The nonphysical zero pseudofermion state |0) gives no contribution to 
the electron scattering, while the two pseudofermion state b^b^\0) can be projected 
out by making it thermodynamically very improbable. This latter can be achieved 
if we choose the pseudofermion chemical potential A to be very high A — > oo. In 
this limit the leading terms to any physical quantity will be dominated by the one 
pseudofermion states and the two-fcrmion contributions will be suppressed by a 
factor ~ e~P x with respect to them. Physical quantities must be normalized to the 
number of pseudofermions present, 2e _/3A . 

The general form of a TLS Hamiltonian can be written as: 

#TLS = A^6+6 Q + J2 A ^ao<V , (2T) 

the Pauli matrices T l aa , (i = x, y, z) denote the 'pseudospin' of the TLS and a = l,r 
labels the states of the TLS. The coefficient A z describes the asymmetry energy 
between the left- and right sides of the TLS while A x and A v stand for the tunneling 
transition (see Fig |l|). The splitting of the lowest-lying two states of the TLS can 
be expressed via these quantities as A = 2(%2 i Af) 1 / 2 . 

As we shall see, only electrons close to the Fermi surface interact effectively with 
the TLS. Therefore, in the spirit of Landau's Fermi-liquid theoryO the conduction 
electrons are treated as noninteracting ones and their Hamiltonian is written as 

H e =^2 ea tn s aens , (2.2) 

where a^ ns and a ens are the electrons' creation and annihilation operators with spin 
s and energy e. The quantum number n classifies the orbital structure of the states. 
It can be, e.g., a spherical wave index-pair n — (I, m) in the free electron case, or a 
crystal field index reflecting the symmetry of the host. 

In our calculations we mostly use a simplified density of states for the electrons, 

= *- D <*< D > (2.3) 

w [0 otherwise, v ' 

where we assumed that the band can be characterized practically by only two 
parameters: the density of states at the Fermi level go and the band width D. The 



role of the energy dependence of the density of states will be discussed in Sec. VII. 

The most general form of the interaction between the electrons and the TLS can 
be written as 

Hi= ]T at ns btV^ a ,b al a^ s , (2.4) 

where /i — 0,2;,?/, z and t®, — 5 aa i. In the following we use the convention that 
Greek indices, /i, v, .., take the values {0, x, y, z}, while the Roman letters, ..,will 



only be used for the components x,y,z. Note that in Eq. (2.4) the couplings 
V^ n , are assumed to be energy- independent. This assumption usually does not 
influence the infra-read behavior of logarithmically divergent models. The real 
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electron spin s plays the silent role of flavor, i.e., the couplings do not depend on it, 
but its multiplicity will essentially affect the ground state and the low temperature 
behavior of the system. 

Assuming a simple two-body interaction b etwe en the TLS and the conduc- 
tion electrons the coupling constants in Eq. (2.4) can be estimated by simple 
integrals .BEia The physical interpretation of these coupling constants is as follows. 
V° ± V z describe the on site electron scattering on the two positions of the TLS. 
If the moving particle is an ion different from the host, then this part of the inter- 
action may be strong, g V° ± QoV z ~ 1. These couplings are responsible for the 
generation of the Friedel-oscillations. The coupling V z , also sometimes called the 
'screening interaction', describes the fact that the electrons are scattered differently 
if the TLS is sitting in its left- and right position, and if the distance d between 
the two sites is small, kpd <C 1, then V z <C V°. It is the coupling V z which is 
responsible for the eventual localization of the TLS discussed in the Introduction. 

The terms proportional to V x ± iV v describe electron assisted tunneling of the 
heavy particle between the two TLS sites, i.e. processes where a conduction electron 
is scattered by the TLS flipping from one side to the other simultaneously. These 
processes are graphically depicted in Fig. |[ They are mostly generated by the 
fluctuations ittfr-thc barrier of the TLS potential due to the local electronic charge 
fluctuations .llao Since they contain the overlap integral between the heavy particle 
states <pi(R) and (f> r (R), they are much smaller then the on site terms: V x ± iV v <C 
V Z M 




(i) 




FIG. 2. Sketch of the screening interaction (i) and the assisted tunneling process (ii). 



III. LOGARITHMIC CORRECTIONS IN LEADING ORDER AND 
ORBITAL KONDO EFFECT 

Having established our microscopic model we turn to the analysis of the pertur- 
bation series. We use the Matsubara-technique to handle the interacting electron- 
pseudofcrmion system. The following imaginary time Green's functions are intro- 
duced for the two types of particlesfl 

G e (e,n,T,e',n',T') = -<T T [a ens (r)aj„, s (r')] ) , 

g(a,T,a',r') = -<T T [6 a (r)6+ (/)] ) . (3.1) 

Here T T is the r-ordering operator, the operators a ens (r) and 6 Q (r) are in Heisen- 
berg representation and the angular brackets denote thermodynamic average. Then 
the unperturbed Green's functions in Fourier representation are given by 

G°( W )=(u,-e)-\ 

g H = (^-AV)-\ (3.2) 
where u> denotes a Matsubara frequency lo = ioj n and the pseudofcrmion frequency 



has been measured from the chemical potential A. In Eq. (3.2) a short matrix 
notation was used for G e and Q whose indices en and a have been suppressed, 
respectively. 
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As in the case of the Kondo effect, the perturbation theories lead to logarithmi- 
cally divergent contributions to different quantities as the relevant energy scale or 
the temperature tends to zero. To demonstrate this we calculate the vertex cor- 
rections of the lowest order depicted in Fig. [| where the solid lines stand for the 
electrons and the dotted lines stand for the pseudofcrmions. Their contribution is 



/ eode ]T (v: n X i n' ^QQi 



f - n(e) 
- e - A*' 



-V J V 1 ,t 1 ( ^ 



(3.3) 



where w is the initial energy of the incoming electron, go denotes the density of 
states at the fermi energy, and n(e) is the Fermi distribution function at tempera- 
ture T. (For the sake of simplicity we put the pseudofermion energy at the chemical 
potential A: lo ps = 0.) 





(i) 



(N) 



FIG. 3. Lowest order vertex corrections in the perturbation theory. 

The energy integral can be approximated by m(Z?/ max{|cj|, T, A}), where A = 

A* 2 ) 1 ^ 2 , and terms of the order of unity have been neglected with respect to 
the logarithm. For the sake of simplicity we will consider the case when |w| is 
the largest energy scale. (The case T > A, \w\ can be treated similarly while the 



discussion of the case A > T, \ ui\ w ill be posponcd to Sec. |V[) Then the vertex in 



the first order approximation (3.3) gets the form 



E*(w) = V? - 2ie ljk g V J V k ln(£/M) + . . . (3.4) 
r°(Lu) = V°, (3.5) 

where e' J is the Levi-Civitta symbol, the underlined quantities are matrices in the 
space of electron orbital states, V^ n , — * V_ l , and summation must be carried out on 
the repeated indices. The logarithmic term occurs only if the commutator V^] _ 
does not vanish (noncommutative model). In the reality, the three matrices V_ x , V_ v 
and V_ z do not commute for a realistic TLS, and a commutative model can be 
relevant only if V_ x and V_ v are negligibly small. We note that the spin i Kondo 
model is also an example of a noncommutative model, where the real spin state of 
the electrons takes over the role of the 'orbital spin' n. 

Continuing the calculation of higher order graphs, higher orders of logarithms 
arise: terms proportional to V n \n n ~ m (D/\u>\), m > 0, in general. According to the 
usual terminology vertex corrections with m = 1 are called 'leading logarithmic' 
while those with m = 2 are referred to as 'next to leading logarithmic' diagrams. 
One can prove that the leading logarithmic contribution is generated only by the 
so-called parquet diagrams - such vertex diagrams which can be cut in two separate 
parts by cutting one electron and^one pseudofermion line in such a way that the two 
parts are also parquet diagrams.t3 In Fig. ||.a we show all the third order leading 
logarithmic diagrams. Up to third order only the diagram in Fig. ||.b gives a next 
to leading logarithmic contribution. Adding the total contribution of the diagrams 
in Fig. ^.a to the ones in Fig. [|we obtain that the vertex function up to third order 
in the leading logarithmic approximation can be written as 



G 



r*(o>) = Yl - 2ie iJ '%y J V fc ln(.D/M) 

- 2^{[Z',F]Z Z + Z'tST.Z^^D/H + ... (3.6) 

At small energies the value of the logarithm can be fairly large \n(D/T) ~ 10 which 
might damage the convergence of the expansion above. The standard way to make 
the sum still somehow convergent is to sum up an infinite number of diagrams in 
the perturbation series that can be achieved by the scaling or renormalization group 
method. The 'Mth order scaling equations' will sum all contributions for which 
m < M. In the present and the following Section we consider the cases M = 1 and 
2, called as leading logarithmic and next to the leading logarithmic approximations, 
respectively. In Section ^ we go even beyond these approximations. 





FIG. 4. Second order vertex corrections generating leading logarithmic (Fig. a.) and 
next to leading logarithmic (Fig. b.) contributions. 

A simple derivation for the leading-4pgarithmic or poor man's scaling was given 
by Anderson for the Kondo problem,a which we now apply for the TLS problem. 
The scaling hypothesis assumes that there exist TLS's having different individual 
parameters V_ l , A% and D but having the same low-energy properties, i.e. they 
have the same effective TLS-conduction electron interaction T(uj). The parameters 
of these equivalent systems are connected by the renormalization group transfor- 
mation. 

The easiest way to construct this transformation is to consider an infinitesimal 
transformation of the form Vj 1 = V 1 + dV_ l , D/D' = 1 + dx, and require the 
invariance of the vertex function ( |3.6| ) for arbitrary values of the dynamical variable 
u. As one can check explicitly, the transformati on D — * D' and V_ l — > V_ l — 
2ie z: ' k goV^V^ In (-§7) really leaves the expression (3.6) invariant. The obtained 



transformation can be cast in the form of the following differential equation: 
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where x = In Dq/D, Do being the initial bandwidth cutoff, and the dimensionlcss 
couplings w M = f?oH M have been introduced. Thes 'leading logarithmic scaling 
equations' have to be solved with the initial condiction that the couplings have their 
bare values at D = Do. The vertex function (|3.6| ) can be obtained by integrating 
these scaling equations from x — to x = In 

One can check that the integration of Eq. ( |3.7| ) really reproduces the higher order 
leading logarithmic terms, and an expliciL summation of the parquet contributions 
leads to the same differential equations!!!! i.e., the scaling hypothesis is correct in 
leading logarithmic order. We remark that the scaling hypothesis is not proved in 
general for its thermodynamic applications, and its validity must usually be justified 
term by term in the perturbation series. 

As we mentioned in Sec. |l[ \v x > v \ <C \v z \ for a realistic TLS. This property 
makes us possible to gain some insight to the low-energy pro perties of a TLS. 



Using the smallness of the assisted tunneling couplings Eq. (3.7) can be linearized 



± iv v to obtain 

dv± 
dx 



= ±2[v ± ,v z ]_ . (3.9) 



and v z is constant in this approximation. Choosing the electron representation 
where v z is diagonal, the solution of the previous equations is given by 

- ^(0)exp[±2x«,„,(0) - < n (0))] . (3.10) 

That matrix element, for which ± (v^, n , (0)— i>* n (0)J is the largest, will soon outgrow 
the others. This gives a two dimensional subspace where the scaled couplings get 
the form of the anizotropic spin | Kondo couplings: 

<n' = wV «n' , (3-11) 

<r l nn , being a Pauli matrix acting in the relevant two dimensional electronic sub- 
space. From this point the problem continues similarly to the anisotropical Kondo 
problemE3 First the smaller couplings v x and v v approach v z , then they grow 
further isotropically and diverge at the energy scale 



< 3 - 12 > 

where v x rand Vq are some characteristic matrix elements of the unrenormalized 
couplings Ji3 Tk is called the Kondo temperature and its index I refers to the 
leading logarithmic approximation. This divergence is clearly nonphysical, since it 
would injply a finite temperature phase transition in an effectively one-dimensional 
systemcil and, as we shall see later, it is only a consequence of the inaccuracy of 
leading logarithmic approximation. 



IV. NEXT TO THE LEADING LOGARITHMIC EXPANSION 



As we have seen in the previous Section the leading logarithmic approximation 
has several discrepancies: It produces an artificial divergency of the effective cou- 
plings at a finite energy scale, and it does not tell anything about the scaling of 
the energy splitting A. As we know from the treatment of the Kondo problem, the 
leading logarithmic results can be considerably improved by going one more loga- 
rithmic order further and both the above-mentioned problems will be solved already 
at the next to leading logarithmic level.Ej Therefore in this Section we extend our 
previous calculation to the next to leading logarithmic order. The summation of 
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the next to leading logarithmic diagrams will be performed in the framework of the 
multiplicative renormalization groupO 

Like in the former Section, the multiplicative renormalization group scheme pro- 
vides connection between equivalent physical systems with different parameters and 
it is formulated as an internal symmetry of the Green's functions. In our case the 
physical systems are characterized by D, and also A 1 is involved because it gets 
perturbative contributions in this order. The main assumption of the multiplica- 
tive renormalization group is that the Green's functions and vertex functions of the 
original and the scaled TLS's have the same functional form, and_they only differ 
by multiplicative factors Z e and Z p , which are independent of 

G e (u;/D', V', A') = Z e (D'/D, V)G e (u J /D, V, A) , (4.1) 
g(u/D', V, A') = Z P (D'/D, V)G(w/D, V, A) , (4.2) 
r{uj/D', V) = Z e (D'/D, VY x Z v [p'lD, Vy^iw/D, V) , (4.3) 

where the primed parameters are those of the renormalized system, and the renor- 
malized electron- and pseudofermion Greens functions are 

G e = (c^-e-EeH)- 1 , (4.4) 
g= (cj-AV-EH) -1 . (4.5) 



In Eqs. (4.1-4.5) the Green's functions and the self-energies are matrices just like 
in Eq. (3.2). In our case one can easily convince himself that in a given diagram 
dressing an internal electron line by a pseudofermion self-energy results in sublead- 
ing corrections in e~^ A . Therefore it follows thatjn the A — > oo limit the electronic 
wave function renormalization factor Z e is one:t2 

Z e = 1 . (4.6) 

As in the p revious Section the infinitesimal renormalization group transforma- 



tions (t4.l| - |4.3| ) can be written as differential equations for the different couplings of 

the scaling equation for the couplings v is 

/?(«) , (4.7) 



the model. The most general form of the scaling equation for the couplings v is 

dv 



dx 

where (3(v) is a polynomial, which is determined by perturbation theory, and the 
different couplings are represented by a single symbol v. In the leading logarithmic 
approximation the polynomial is of s econd order, and it can be determine d by 



taking the derivative of the vertex (3.2) with respect to In D. In this case Eq. (4.7) 
generates the leading logarithmic terms to all orders. To get the next, third order 
term in (3{v) the vertex and pseudofermion self-energy corrections must be treated 
simultaneously. 




The get some insight what the physical meaning of the scaled couplings is, let 
us assume for a moment that the temperature is the largest energy variable of 
the vertex function (i.e., we investigate scattering of electrons at the Fermi surface 
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off the impurity), and apply the renormalization group transformation with the 
special choice D' — T. Then, since In Sr m no vertex corrections appear for the 

scaled TLS with the primed parameters. Therefore r*(T, D',V) « V'\ Thus the 
solution of Eq. (4.7) at the scale x = ln-p, up to the numerical factor Z p , is just 



the effective interaction of the TL S an d the conduction electrons at temperature 
T. (Knowing the solution of Eq. (17) one can easily determine the factor Z p as 
well). Similar procedures can be applied to calculate -fthjffiical quantities like the 
electronic scattering rate or the impurity specific heat.E3o 

Calculating the vertex and pseudofermion self-energy diagrams shown in Fig. ||, 
we obtain 



f = v l + {-2is ijk (v j v k ) + N f [2v 3 Tr(uV) - v l Tr(uV)]} H D /M) , ( 4 - 8 ) 
E(w) = -NfTrivVyluiI- AV=]t j ' ln^D/M) , (4.9) 



where 7 = QqT is the dimensionless vertex, the trace operator Tr{...} is acting in 
the electronic indices, and Nf is the spin degeneracy of the conduction electrons 
(Nf = 2 for the physical case). We remind the reader that underlined quantities 
denote matrices in the orbital indices of the electrons. The self-energy contains 
terms proportional to win — , which give a contribution to Z p . 

To construct the scaling equations, an infinitesimal step D/D' — 1 + dx 
(x = MDq/D)) must be considered in Eqs. gQ, (|J), Q and Q. Then 
plugging Eq. (4.9) into Eq. ([D]) the renormalization factor Z p can easily be deter- 
mined by comparing the ~ 1/ui terms on the two sides of Eq. (4.2) : 



Z p (D'/D, V') = l- N f \n(D'/D) Tr^V ) 



(4.10) 



Kno wing Z p o ne can also read off the infinitesimal transformations of A 1 from 
Eq. Q and (jl|): 

(4.11) 



A' 1 = A 1 + 2N f ln(D'/D) [A* Ti{v!v j ) - A j Tr(^w J ')] + 0(v 3 ) 
which can be written in a differential form as 
dA i 



dx 



-2N f [A* TrfeV) - A 3 Tr(wV 



(4.12) 



This equation describes the renormalization of the TLS splitting as a function of 
the temperature. One can easily prove that the A l 's are always decreased under 
the scaling transformation. 

The scali ng eq uations for the couplings v l can be generated in the same way 
using Eqs. fl4T0|) , (|J), and (|J): 



— = -2ie ijk v j v k - 2N f \v l Tr(vV) - v j Tr(v j v i 
dx ~ 



(4.13) 



These scaling equations can be solved with some realistic restrictions on the initial 
couplings and it can be shown just like in the leading logarithmic case that the 
relevant-electron subspace is two-dimensional in the weak coupling regime of the 
scaling.E3 As we will show in Sec. ^ t his st atement remains valid even below Tk 
and the only stable fixed point of Eq. ( [4.13 ) is where the interaction has the same 
structure as the spin \ Kondo coupling, v l nn , ~ a l nn , . 
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FIG. 6. Scaling of the dimensionless couplings and the splittings for a simplified TLS 
with v l — v l a_ l in the leading- (dashed lines) and next to leading logarithmic (continuous 
lines) approximation. 

Typical scaling curves are given in Fig. ^, for the special case where the u"s 
are proportional to some Pauli matrix w l = nV ! (no summation). As it is shown 
in Ref. |l8| this is a fairly good approximation for a TLS interacting with a free 
electron band. Both the leading logarithmic and the next to leading logarithmic 
scaling trajectories are presented. As one can see from the Figure, due to the 
presence of the third order terms in the scaling equations the couplings v l do not 
diverge anymore, but they scale to a finite value. 

Then the orbital Kondo temperature is identified by the characteristic energy 
scale, where the couplings approach their fixed point value, and it can be expressed 
as 



'°('«>*" 4 (4r . 

where the indices II and refer to the next to leading logarithmic approximation 
and the initial value of the parameters, respectively. This Kondo scale like is 
invariant under scaling. 

Due to the appearance of the prefactor this Kondo temperature is considerably 



smaller than the one obtained in the leading logarithmic approximation, Eq. (3.12). 
While for reasonable parameters T£ is of the order of ~ \0K in the next to leading 
logarithmic approximation TU ~ IK — O.IK. 

To close this Section we estimate the renormalization of the tunneling amplitude 
A x . To be specific we consider a symmetrical TLS with A z = A y — 0, and the 
above-mentioned simplified couplings v l — uV\ Then the scaling equation for A x 
can be written as 

l^!=-4AT / [(^ + K) 2 ] . (4.15) 

The scaling procedure should be stopped by A x for very small x, where the splitting 
becomes the dominating l ow- energy scale in the argument of the logarithms (see 



the discussion below Eq. (3.3)). The corresponding energy scale is defined by the 
implicit equation A x (ln(Do/T*)) — T*. Below this energy scale the motion of the 
TLS is frozen out and remains constant for D < T*. Taking into account 



that in most of the scaling procedure v z ss Vq ^> v x in Eq. (4.15), the 'freezing 
temperature' T* can be roughly estimated as 

T* = A -M . (4.16) 



D 
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T* is-sometimes also referred to as the renormalized or effective splitting of the 
TLS.E3 For realistic parameters T*/Ag can be as small as 10 2 . Taking into ac- 
count, however, that in the noncommutative model v z and v x also depend on x one 
obtains that A x can be reduced by several orders of magnitude for large enough v z 's 
as can be seen in Fig. ^|. A more detailed discussion shows that the renormalization 
of A z is usually much smaller, and the more symmetric the potential the larger is 
the reduction of energy splitting. 



V. STABILITY OF THE TWO-CHANNEL KONDO FIXED POINT: 1/N F 

EXPANSION 



As we have seen in the previous Sections the noncommutative TLS shows a 
strongly correlated behavior at low temperature. This behavior manifests in the 
appearance of logarithmic singularities in the vertex function which can be associ- 
ated to the formation of a Kondo-like ground state with a binding energy of the 
order of Tk- Using the multiplicative renormalization group method we were able 
to handle these logarithmic singularities. In the next to leading logarithmic ap- 
proximation the artificial divergence of the vertex function has been removed and 
we could also account for the renormalization of the TLS splittings. However, for 
small values of Nf (including the physical case Nf = 2) the next to leading loga- 
rithmic approximation breaks down i n the vicinity of the fixed point, as the higher 
order terms of the j3 function in Eq. (4.7) become comparable with the first ones. 
The main purpose of the present Section is to go below the Kondo scale, and see 
what happens to the TLS at energy scales T,w < Tk- We shall circumvent the 
above-mentioned difficulty by making a 1/Nf expansion. 

We have seen in the previous Sections that at high temperatures there are two 
orbital electron channels which dominate the scattering. The logarithmic anomalies 
are essentially due to these two channels and the other orbital channels give only 
a small contribution. We have also shown that if only these two channels are 
considered, then in the weak coupling region the (now 2x2) matrices v z scale 
towards simple spin 1/2 operators v % ~ a 1 = 2S\ with S e = l/2.olia Gided by 
these facts it has been conjectured in Ref . [l8] that at the low energy fixed point the 
electron-TLS interaction is described by a simple effective exchange interaction in 
the orbital degrees of freedom ~ t 1 SI- and thus the model is equivalent to that of 
the two-channel spin Kondo problemu which exhibits non- Fermi liquid properties. 
The double degeneracy of the two 'channels' in the latter model corresponds to the 
real spin indices of the conduction electrons in the TLS case, and the splitting A 
of the TLS acts as a local magnetic field at the imputity site in the two-channel 
Kondo analogy-iSfi£-also Table I an Sec. ^). Several recent experiments on metallic 
point contactsOC3o have been interpreted in terms of the complete equivalence of 
the aforementioned two models (See also Sec. VIII). This 'complete equivalence' 
is, however, far from being trivial. As indicated above, the mentioned results of 
Vladar and Zawadowski are only trustworthy in the high-energy (weak coupling) 
region where the renormalized couplings are small, and their approximations loose 
their meaning in the strong-coupling region T, u> ~ Tk ■ 

Vladar and Zawadowski have already remarked that the TLS model has other 
fixed points which differ from the previously mentioned one by the number of orbital 
channels coupled to the TLS. An example of such fixed points is given by v 1 ~ 
SI with S e > 1/2. Simple estimations show that the coupling of the different 
orbital electron channels is strong thus one would naively expect that as soon as 
the couplings of the two orbital channels dominant at T ^S> Tk become of the order 
of unity all the other couplings start growing up as well. Therefore, considering 
for instance three different orbital channels one could imagine that although there 
are two dominant channels at high temperatures (T > Tk), the third channel 
becomes also important below the Kondo temperature Tk-, and the low temperature 
scattering is described by an S e — 1 orbital electron spin corresponding to the 2S e + 
1 = 3 orbital channels. Thus, naively, one could expect a series of Kondo effects 
corresponding to the increase of the couplings of the different orbital channels. 
Vladar and Zawadowski have also remarked already in their early worklla that 
below Tk higher angular momentum scattering might be relevant. 
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Furthermore, as it is wellknown from renormalization group theory^! the low 
temperature behavior of a model is not only determined by the structure of the 
stable fixed points it scales to, but also by their operator content. Therefore, to 
identify two models one must be very careful and has to consider the operator 
content of the two models as well. There are several models, which possess a 
simple fixed point, but have nontrivial operator content. An example is the fixed 
point found by H. Pang for a generalized two channel Kondo problcmcJ which 
has a spectrum composed from two independent Fermi liquid spectra, however, its 
thermodynamical properties are claimed to be determined by a non-Fermi-liquid 
leading irrelevant operator. 

Our purpose in this Section is to establish a more rigorous correspondence be- 
tween the two models and to show that the stable low-energy fixed point of the 
model is correctly described by the conduction electron orbit al s pin S e = 1 /2. For 
this purpose we investigate an A^-flavor TLS Hamiltonian (2.2) and (2.4) with a 
general spin degeneracy Nf and develop a systematic l/Nf expansion. Then we 
classify all the possible fixed points of the model and investigate their stability and 
operator content. We find, that the only stable fixed point is the S e — 1/2 fixed 
point independently of Nf and the number of orbital channels considered. While 
we find that the operator content of the A^-flavor TLS model is much richer than 
that of the A^-channel spin 1/2 Kondo model, the dimension of the leading irrele- 
vant operators, and thus the thermodynamical and dynamical behavior of the two 
models is qualitatively the same. Since our result are exact in the limit where Nf 
is large and they are independent of the special value of Nt we expect that these 
results are also valid even for the Nf = 2 case. 

In these considerations we first neglect the role of the splitting of the TLS which 
serves as a lower cotoff in the scaling procedure. The role of this splitting will 
be discussed later. First we describe the calculation in the l/Nf order. Then we 
discuss the technicalities of the extension of the calculation to l/Nf order, and 
finally we exploit the mapping found to describe the non-Fermi liquid properties of 
the TLS in terms of the 2-channel Kondo model. 



A. The l/Nf order analysis of the scaling equations 

The next to leading logarithmic scaling equations for the A^-flavor TLS model 



have been derived in Sec. IV 
dv 



/ = -2ie ijk v j v k - 2N f v l Tr {v j v j } + 2N f v j Tr {uV } . (5.1) 



We remind the reader that the Nf factor appearing in the last two terms of Eq. (5.1) 
are due to the presence of the electron loops in the third order vertex correction 
and the pseudofermion self energy in Fig. a. We stress that the scaling equation 



(5.1) is very general. To derive it one has to assume only that some sort of Fermi 
surface exists and that the conduction electrons can be described in terms of the 
Fermi liquid theory. The special choice of the density of states has no effect on the 
universal properties of the sys tem up to this order. 

The fixed points of Eq. ( p.l| ) are determined by the condition that its right-hand 
side vanishes. From the presence of the factor Nf immediately follows that the 
fixed point couplins of the TLS should be of the order of ~ l/Nf. Therefore, as 
observed first by Nozieres and Blandin for the overscreened Kondo modcljj in the 
Nf — > oo limit the scaled couplings remain in the small coupling region and the 
scaling equations give exact results. (This should be contrasted to the Nf = 1 
'simple Kondo' gase where it can be shown with other methods that the couplings 
scale to infinity.Q) 

The fact that the fixed point coupling is of the order of l/Nf makes alsapceaible 
to develop a systematic l/Nf expansion for our A^-flavor TLS model. c3'El'E3 To 
understand this point let us investigate an n'th order logarithmic vertex correction 
in the perturbation series. Around the fixed point each vertex is of the order l/Nf 
which results in a factor l/Nf. This small factor is partially compensated by the 
electron loops present in the diagram bringing up a factor Nf , where L denotes the 
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number of electron loops in the diagram. Since this last number is always smaller 
than n/2 up to a given order in 1/iVy there exist only a finite number of diagrams 
which may generate corrections to the scaling equations around the fixed point. 
This analysis can be extended to terms generated by the pseudofermion self-energy 
corrections. Therefore, if one wants to calculate the fixed point coupling or the 
scaling exponent of the different operators up to a given ord er in 1 /Nf only a finite 



number of diagrams must be considered (see also Subsection VB). A quick analysis 
shows that the lowest order diagrams in 1/Nf are exactly the ones generating the 
next to leading logarithmic scaling equations shown in Fig. |j| 



Now we turn to the analysis of the classifi cati on of the fixed points of Eq. (5.1). 
It is easy to show that the last term in Eq.( |5.ip can be eliminated from the fixed 
point equation by making an orthogonal transformation v l — > Y) , OijV? , O being 
an orthogonal matrix. T heref ore, it is enough to consider the first two terms on 
the right-hand side of Eq.(|5.l|): 



e ijk v j v k = iN f v l Tr {v j v J } . (5.2) 



Multiplying Eq.(|5.2|) by v l and taking its trace one obtains the following equations: 



where cm = Tr{v l v 1 } (i = x,y,z) and f3 = Tr{v x v v U z — v z v y v x }. From Eq.( |5. 
immediately follows that either at least two of the on 's are zero at the fixed point 
or they are all equal: a x = a y = a z = a. The first case corresponds to the 
commutative TLS, where the assisted tunneling is ignored and the couplings v_ x 
and v y identically vanish. This is evidently an unstable fixed poi nt. In the second 



case it is worth introducing the matrices J 1 = 2 N fa — ■ Then Eq.(5.2) tells us that 
the J"s satisfy the SU(2) Lie algebra 

[J\ J j ] — ie ljk J k . (5.4) 

Therefore the general form of the J l 's at the fixed point can be given by a direct 
sum of finite dimensional irreducible representations of the SU(2) spin algebra 

n 

J* = @ Sfy , (5.5) 

where the 's denote integer or half-integer spin representations and n is the 
number of irreducible representations involved. The S^'s in Eq.( |5.5| ) are acting 
only in a finite dimensional subspace of the total electronic phase space, and in the 
rest of the phase space the J l 's give identically zero. The value of a can easily be 
determined using its definition: 

(5.6) 



^NjY2 k S ik) (S {k) + l)(2S ik) + l) 



It is obvious from Eq.(5.6) that the value of a depends in an essential way on the 
electronic orbital spin structure of the fixed point approached by the w"s. Carrying 
out the scaling all the a,'s become equal and they scale to a special value of a 
characteristic to the fixed point. 

The real low-energy properties of a TLS can only be described by some stable 
fixed points. The presence of the other unstable fixed points could only be observed 
if the paramaters of the TLS are finetuned- but they cannot produce a universal 
scaling of the diff erent physical quantities.cH Therefore, having found all the fixed 



points of Eq. (54), we turn to the stability analysis of the fixed points. 



1. The unstable fixed points 

We first show that the only possible stable fixed point is the one which is equiv- 
alent to a single S e = 1/2 orbital electron spin. The proof proceeds in two steps. 
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First we show that the 'composite' fixed points where the representation J 1 is re- 
ducible and contains several different spin representations are unstable. In the 
second step we prove that the fixed points where the J l 's are irreducible but are 
equivalent to an S e > 1/2 spin representation are also unstable. Therefore, one can 
conclude that only that fixed point can be stable where the couplings yf contain a 
single S e = 1/2 spin representation. The stability of this latter will be proved in 
the next Subsection. 

To prove the instability of a fixed point our strategy is to find a single unstable 
scaling trajectory running out from the fixed point. If the v l, s are composed from 
at least two spin representations we can chose two of them, 5m and 5( 2 ), and 
consider small deviations from the fixed point vf — (2NfCt) J2 k ^Yfc) °f the f° rm 

6v* = 2N f (6a 1 S i il) + <fa 2 5[ 2) ) . (5.7) 

These special scaling trajectories lie in the subspace expanded by the two arbitrarily 
chosen spin representations Sm and 5(2). Linearizing the scaling equations one 
obtains a closed system of equations for the 5a>i 's which can be solved easily. Then 
it is trivial to show that the operator 



(aiSfo - a 2 S\ 2) ) (5.8) 



with 1/aik = ^NjS [k) {S [k) + l)(25 (fe) + l)/3 scales like - e 4N ' ax ~ D~ m f a and 
is relevant at small energy scales (low temperatures) D — > 0. Therefore, any 'com- 
posite' fixed point is unstable. 

Now we proceed by proving that the 'irreducible fixed point' v l ~ 5 l is also 
unstable for 5 > 1. To prove this it is enough to find a relevant operator in the 
space of the Sv_ l, s. A detailed discussion of the construction of such an operator is 
given in Ref. plj here we only give its explicit form: 

0rel ~ 6a(2N f ){-lwS- + 5*5*) + r^S* + S»S')] 

+ t z (S x2 + S y2 -2S z2 )} . (5.9) 

Substituting this expression into the scaling equations we find that O tc \ scales like 
~ T- SN f a and is relevant. 



It is obvious from Eq.(5.9) that the above considerations break down fo r an 
5=1/2 orbital electron spin, since for 5=1/2 the operators appearing in (|5.9|) 
vanish identically as a consequence of the special properties of the 5 = 1/2 spin 



algebra. Thus, we can conclude that any low temperature fixed point of Eq. fl5.lD 
which is composed from several spins or corresponds to an 5 > 1 spin representation 
is unstable. 

2. Stability analysis of the 5 e = 1/2 fixed point 

We now prove that this fixed point is really stable. To examine the stability we 
chose a basis in which the i/'s at the fixed point can be written in the blockmatrix 
form 

{v - ] * = W] ( o o) ' ( 5 - 10 ) 

where a 1 (i = x,y,z) denote the Pauli matrices. For the sake of simplicity we 
assume that only a finite but arbitrarily large number N Q of orbital channels are 



considered, thus the v l, s appearing in Eq.(5.10) are N a x N a matrices. Then the 



deviation of the v l, s from their fixed point values can be written in the form 

**-(££)■ (5 ' 11 > 

where rf, M_ l , f and t i+ denote 2x2, (N Q - 2) x (N Q - 2), 2 x (N a - 2) and 
(N — 2) x 2 matrices, respectively. The matrices rf and M' are Hermitian while 
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t l and t l+ are Hermitian conjugates of each-other. The linearization of the scal- 
ing equations is lengthy but straightforward and one obtains the following set of 
equations: 



dM l 
dx 
drf 
dx 



dt l 
dx 



2 

-4r Ml 



(5.12) 



j,k j^i 

E^^F7 Tr W+^V}, (5.13) 

(5.14) 



__L^ e «vt*-— t 

N f ^ ~- Nf~ 

J j,k 1 



is rather 



The matrices f satisfy the Hermitian conjugate of Eq.(5.14). 

The detailed stability analysis of the linearized scali ng e qu ation s 
tedious.cil However, due to the complete decoupling of Eqs.( pTl2|) f ( 5.12 ), and (5.14) 
all the scaling exponents can be calculated exactly. Eqs. Q5.12] — 5.14) have an 



infinite number of marginal solutions remaining unsealed under the linear approxi- 
mation. A thorough analysis reveals that these marginal operators do not affect the 
SU (2) spin algebra of the matrices (v^fp and they correspond to the rotations of the 
pseudospin of the TLS or the two-dimensional electronic orbital subspace, where the 
SU(2) algebra is realized.^ All the other operators scale like ~ e~ Ax ~ (D/T K ) A 
to zero as D — > with an exponent A > and are irrelevant at low energies. 




FIG. 7. Sketch of the attractive S e = 1/2 manifold embedded into the space of the 
couplings v_ l (i = 1, 2, 3). The S B = 1/2 fixed manifold is represented by a heavy contin- 
uous line while the scaling trajectories are indicated by arrows. A marginal perturbation 
corresponds to moving along the heavy continous line. 

Accordingly, the S e — 1/2 'fixed point' of the scaling equation is rather an 
attractive 'fixed manifold' embedded in the manifold of the general couplings v l . 
This situation is sketched in Fig. |7[ The unitary transformations of the S e = 1/2 
spin algebra are associated to trajectories lying in this attractive fixed manifold and 
correspond to zero exponents. The number of leading irrelevant operators is also 
infinite. These leading irrelevant operators associated to the exponent A = 2/Nf 
are given by the following expressions: 



i 



(5.15) 
(5.16) 



Or 



(5.17) 



where the Q l, s and J l 's denote arbitrary small constants, and a is the unit matrix 
acting in the 2x2 block of the matrices in Eq.(5.10). The operator Om appears due 
to the presence of the other orbital channels of the conduction electrons. The second 
operator, Oj corresponds to the leading irrelevant operator of an anisotropic Nf 
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channel Kondo problem .eIHHSeS Finally, the third operator describes the splitting 
of the TLS generated by the two dominant electron channels coupled to the heavy 
particle. 




0.0 10.0 20.0 30.0 40.0 50.0 0.0 10.0 20.0 30.0 40.0 50.0 



ln(D /T) ln(D„/T) 

Fig. |a Fig. |b 

FIG. 8. The scaling trajectory of the matrix norms at = Tr{i/ } (i = x,y,z) for the 
Nf — 3 case (Fig. and their relaxation to their fixed point value (Fig. ^|b) for Nf = 2 
(curve a.), Nf = 3 (curve b.) and Nf — 4 (curve c.) in a logarithmic scale. 

The statements above have also been tested numerically. In Fig. || we show the 
scaling of the norm cti = Tr{(tf) 2 }. The initial values of the couplings have been 
estimated by assumingj-ar-screened Coulomb interaction between the conduction 
electrons and the TLSii§El As one can see these scale to their fixed po int value 
a = with the power law dependence ~ D 2 l N s in agreement with Eqs. ( |5.6| ) and 

the scaling dimension of the leading irrelevant operators. The Kondo effect can be 
identified as the breakdown of the curves. 

B. Scaling analysis up to the order 1/Np 

To obtain the 1/-/VJ order scaling equations one has to determine the vertex func- 
tion and the pseudofermion self-energy up to the orders l/iV§ and 1/Nj, respec- 
tively. The corresponding self-energy and vertex diagrams are shown in Figs. ^| and 
|To| . They can, in principle, be calculated easily using Abrikosov's pseudofermion 
technique, however, there are some technicalities which are crucial to obtain the 
correct results. 




e- f. 
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FIG. 9. The pseudofermion self-energy corrections up to the order ~ 1/Nt. 
and the continuous lines denote the pseudofermion and the conduction electron propaga- 
tors, respectively. The cross indicates the contribution of the counterterm, which must be 
calculated up to the order ~ 1/iV/ . 

The first poblem one has to face is the appearance of spurious divergencies in the 
perturbation series. The second order self-energy diagram, Fig. ||a contains a log- 
arithmic contribution but it also gives a constant term <5A = — 2 D In 2 NfTr{v_ l v 1 } 
which renormalizes the chemical potential ('mass') of the pseudofermions. One 
knows from quantum field theory that such terms proportional to D are artificial 
and can be scaled out by introducing (up to second order) a counterterm in the 
Hamiltonian 

tfcount = -5\J2b+b a . (5.18) 

a 

This counterterm is then canceling the spurious divergencies in the different dia- 
grams, and guarantees that the pole of the pseudofermion Green's function remains 
unshiftcd. It has to be determined order by order up to the desired order accuracy 
in the perturbation theory. The inner loop in the self-energy diagram of Fig. ^|.d, 
e.g., contains a constant part which results in an artificial contribution to this dia- 
gram. This artificial contribution is then canceled by the diagram |[e. If one uses 
this counterterm procedure then the Abrikosov projection^ of the b+b a = 1 
subspace should be carried out with the original A parameter in the Hamiltonian. 

Another crucial problem is the separation of the different logarithmic terms in the 
perturbation series. In the perturbative expansion each term of the perturbation 
series can be expanded in terms of the polinomials u> n ln m (u> / D) where ui represents 
the relevant energy variable. In a similar, way the total vertex and self-energy 
functions can also be expanded in terms of such polinomials, and they can be 
written schematically as 

71,771 

r(w) = T nm u n ln m (oj/D) , (5.20) 
where the coefficients are some complicated matrix functions of t he c ouplings. Put 



in another way, the renormalization group hypothesis Eqs. (4.1-4.3) assumes that 
one can change D —> D' and y_ l —>■ v n in such a way that the expansion above 
remains invariant up to a constant multiplicative factor. This holds to all the terms 
appearing in the expansion, which are all, we stress again, well-defined functions of 
the couplings. If one consideres the r 2 hi 2 (w/.D) vertex corrections, e.g., rescaling 
D one generates terms of the type ~ \n(uj/D) ln(D/D'), that renormalize only 
coefficients like Toi. The scaling equations can be, of course, obtained equally 
well from the comparison of the terms Toi and To2- These latters are, however, 
much more difficult to calculate than the coefficients Too and Toi. Therefore the 
easiest way to generate the scaling equations is to collect the ~ ln(D/u>) vertex 
corrections only. The difficulty appearing now is to separate such contributions. 
The contribution of the self energy diagram Fig. ^|.d together with the regularizing 
counterterm contribution ^.e, e.g., can be calculated as 

Egd+e = -uj (~ In 2 -R- + QB.2-2) In -^-\ N 2 f (Tr{uV}) 2 + cnst. , (5.21) 
and it contains both lo In 2 u> and to In u> contributions from which only the second 



one contributes to the wave function renormalization factor Z p in Eq. (4.2) and thus 
the scaling equations. This example shows, that in order to get the right scaling 
equations one has to determine all the subleading contributions very carefully. It 



also demonstrates, that Eqs. (4.1-OJ) and ( 5. 20] ) must be taken very seriously. 
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FIG. 10. The vertex corrections up to the order ~ 1/-^/- F° r the sake of simplicity 
only diagrams without counterterm correction are shown. The missing diagrams can be 
generated by reversing the pseudofermion and the electron lines. 

Apart from these subtilities the calculation is straightforward and after a rather 
tedious calculation the inverse pseudofermion Green's function can be expressed as 
(T = 0): 

g- 1 = uj{1 + (1 - ln2)N f O" + ln(D/u)[N f O jj + 12N f fi 

- (5 - 3ln2)N]0"O kk - (4 In 2 - 6)NfO jk O kj ] } , (5.22) 

where O 1 ^ — Tr(v_ l y?), (3 = —i Tr(v x v y v z — v z y?v x ) and a summation must be car- 
ried out over the repeate d indi ces. Only the leading term proportional to ln(Z)/cj) 



is given explicitly in Eq. ( 5.22 ) and we have taken into account that the couplings 
are of the order of ~ 1/Nf. The leading terms of the vertex function are given by 

r = v l - \n.2Nf(20 ij v j - C^'V) - ]n(D/u)[2ie ijk v j v k - N f (20 lj v> - J V) 

- 4(2 - ln2)N f ie jkl O lj v k v l + 2\n2N } ie ljk O u v j v k + (2 + 5 In 2)N 2 f O"O kk v i 

- (8 + l2]n2)N]O kk O ij v j + (8 + 12 In 2)NfO ij (y k v k 

- (2 + Aln2)NjO jk O kj v!] , (5.23) 

where e ljk denotes the Levi-Civita symbol. 

Having determined the vertex function and the pseudofermion Gre en's fun ction 
we can generate the scaling equations by plugging them into Eqs. ( |4.1| - |4.2| ) in a 
selfconsistent way and we finally obtain 

tin' 

= 2ie ijk v j v k + 2N f (0 3j v l - O lj v°) - 8(1 - ln2)N f ie jkl O i: >v k y! 

+ 16 In 2 Nffiv! + 8N]O ij O jk v k - ( )'■' r' , (5.24) 

where the scaling variable x — ln(D/D') has been used. One can easily check that 
with the substitution v l — > ja l /4 these equations simplify to the scaling equation 
derived by Gan et al. for the iV/-channel Kondo modelo 

The stability analysis follows exac tly t he same lines as in the previous Section. 



The only stable fixed poipt of Eq. ( 5.24 ) is identical with the fixed point of the 
iV/-channel Kondo model:Eil 
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where the Pauli matrices <7 l are acting in a two-dimensional subspace in the elec- 
tronic orbital space. In order to determine the leading irr eleva nt operators one has 



to linearize the scaling equations around the fixed point (5.25). For this purp ose it 
is convenient to write the deviation from the fixed point again in the form ( 5.11 ) 
and then the linearized equations become: 



dM l 
dx 

d£ _ _ i_ 
~dx" 2 



-AM' , (5.26) 
A £ ijk {v j e k + l 3 Z- k ) - A i z (5.27) 

^=--AV e «V* fc - Af , (5.28) 
ax 2 

where the critical exponent A = (2 — A/Nf)/Nf is characte ristic to t he stable 



fixed point. These equations have the same structure as Eqs. ( 5.12 - 5.14 ) the only 
difference is the appearance of the exponent A instead of the factor 2/Nf. Therefore 
all the conclusions concerning the operator content of the fixed point remain in 
princ iple unchanged. It is worth mentioning that, while the fixed point couplings 



( 5.25 ) contain a factor In 2 reminiscent of the special band structure used, there is 
no similar factor in the expression of the critical exponent A, which is universal. 

Similarly to the l/Nf -order case, there exist an infinite number of leading ir- 
relevant operators, in contrast to the multichannel Kondo model. Theoretically, 
these new leading irrelevant operators could change the low temperature properties 
of the model compared to the Kondo model. However, calculations of different 
measurable quantities indicate that while they give a contribution to the physical 
quantities, they do not change their critical exponent. For the free energy, for ex- 
ample, we find a ~ j^a+i behavior, while the scattering rate of the conduction 
electrons scales as ~ T A . With elementary considerations one can also derive the 
scaling of the 'impurity magnetization', — dF imp /dA and the 'impurity suscep- 
tibility' xa = d 2 F imp /dA 2 as a function of the splitting A and the temperature T. 
All these results are in agreement with the eptact ones obtained for the iV/-channel 
Kondo model discussed in the next SectionO'Lj However, it would be interesting 
to find some measurable quantity (like the Wilson ratio,e3 for example) where the 
presence of these operators is manifest. 

Naturally, the l/Nf expansion breaks down in the physical limit, Nf — > 2. How- 
ever, the structures of the leading irrelevant operators being independent of the 
value of Nf we think that — sim il arly t o the i Vf-ch annel Kondo modeled — they 



are correctly given by Eqs. (5.15), (5.16), and (5.17). Therefore we think that the 
properties of a physical TLS with Nf = 2 could be determined by conformal field 
theory methodsEa using a two-channel Kondo model, where these leading operators 
are also included as perturbations around the fixed point. We also expect that the 
scaling exponent of these operators remains degenerate even for Nf — ► 2 and . 
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for Nf = 2 they coincide with the exact Kondo exponent, A exact = 2/ [Nf + 2). I 



C. Analogy with the 2-channel Kondo model and scaling behavior 

In the previous Sections we have basically established a mapping of the TLS 
model to the iVy-channel Kondo model. We have show n that at the fixed point 
the couplings v l take the particulary simple form ( 5.25| ), v l nn , = va % nn , : where the 



indices n and n' label the two relevant orbital (angular momentum) channels. Thus 
the TLS effective Hamiltonian can be written as 

H int = « E a ^^C'C' 6 «' fl eV, ■ (5.29) 

«, a, a> r 
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As the calculation in Sec. IV demonstrates, the TLS Hamiltonian becomes isotrop- 



ical in the couplings v x , v_ v , and v z even above Tk- Moreover, we have shown in the 
previous Subsections that the leadi ng irr elevant operators of the TLS have the same 
scaling dimen sion as the operator ( 5.29 ). Therefore we conclude that the effective 
Hamiltonian (5.29) is adequate to describe the TLS at temperatures around and 
below Tk- 



Now the main observation is that the Haxniltonian (5.29) is formally equivalent 
to the one of the A^-channel Kondo modeE 



H 



Kondo 



N, 



f— 1 e,e' ,<T,cr f 



o-,cr' a ecr-J U, £'<?';f 



(5.30) 



where now the 1/2 inpurity spin is coupled to the real magnetic spins {a — ±) of 
Nf independent electron channels via an antiferromagnetic exchange interaction. 
The correspondence between the two models is explained in Table I. 



multichannel Kondo model 



TLS model 



Hi = JY, sl <,cT> a tvj a e'<T';f H i tt = v Y, a tn S bi<J l nn ,T l aa ,b al a e > n r 



impurity spin S" 1 

electron spin a 
flavor /= l,..,N f 
local magnetic field h 



TLS pseudospin r* 
electrons' angular momentum n = s,p, 
electrons' real spin s = 1, .., Nf 
TLS splitting A 



TABLE I. Correspondence between the multichannel Kondo and the TLS model. 



The multichannel Kondo model has aJ^ ea | dy ,-b eeB ^ fcj ^i ti d .by a variety of methods 
and its properties are well-understood.B^K'y'cJ'Ea'EaOEa Therefore our analogy 
between the two models is particularily useful to understand the behavior of a 
TLS-model in details and translate the exact results obtained for the multichannel 
Kondo model to the TLS case. 

For the sake of simplicity let us discuss first the case, Nf = 2. As we discussed 
above the physical TLS model is mapped to the two-channel Kondo model, which, 
in contrast to the usual single channel Fermi liquid Kondo model, belongs to the 
class of overcompensated spin models and exhibits a non-Fermi liquid behavior. The 
origin of the non-Fermi liquid behavior can easily be understood as follows. In the 
single channel CclSC, clS the temperature is lowered below Tk, a conduction electron 
is bound to the impurity spin by the antiferromagnetic interaction, and they form 
together a spin singlet state (see Fig. f[l].a). This composit singlet without inner 
degree of freedom is a 'dead body' for the rest of the electrons, which are only 
slightly scattered by it, and therefore the ground state of the system is properly 
described by a Fermi liquid theory. 

On the other hand in the two-channel Kondo model around Tk both channels cou- 
ple to the impurity spin antiferromagnetically and therefore, as shown in Fig. [TT|.b, 
they form a bound state with an effective spin 1/2 pointing to the opposite direc- 
tion as the impurity spin. A careful analysis shows that the resulting bound object 
is coupled again antiferromagHetically to the rest of the conduction electrons, and 
generates a new Kondo effect. a Thus, in this picture the non-Fermi liquid behavior 
can be viewed as a never-ending series of Kondo efects at different energy scales. 
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Fig. b. 



FIG. 11. Comparison of the single channel Kondo model and the two-channel Kondo 
model. In the two-channel Kondo model case the impurity spin is overcompensated by 
the conduction electrons. 

The fact that the conduction electrons cannot screen the impurity spin in the 
multichannel Kondo model is reflected in the appi^afl«M>6$ of a nonzero residual 
entropy due to the degeneracy of the ground state:! 



AS, 



(2CKM) 



fc B -ln2 
2 



(5.31) 



Eq. ( 5.31 ) also gives the correct entropy for the TLS .Model in the temperature 
range T* < T < Tk, T* being the freezing temperature.c3 However, for a real TLS 
the splitting acts like a local magnetic field in the two-channel Kondo model, and 
splits up the degeneracy of the ground state, therefore the splitting drives finally 
the TLS towards a Fermi liquid state and the residual entropy vanishes due to the 
presence of the splitting as T — > 0: 



ACp S) = 



(5.32) 



From the theory of the two-channel Kondo model we know that .. th e specific 
heat and the different susceptibilities show logarithmic anomaliescaElJO and that 
the resistivity scales like T 1 / 2 below Tr-O These results can be translated for the 
TLS case without difficulty keeping in mind the analogy (Table I.) of the different 
quantities. The behavior of the corresponding physical quantities for the TLS model 
is summerized in Table II. both for the cases Nf — 2 and Nf > 2, where the relevant 
energy variables, A and T, are assumed to be larger than the freezing temperature 
T* but smaller than T K . 



Nf>2 



^ 2 v ' 

(2-N f )/(N f +2) 



Nf 



xi p (A,T = 0) 
xi p (A = 0,T) 



--imp 



Afli, 



IT 



(*) 



(*) 



(2-^)7(^+2) 



2/(2+N s ) 



In (A/TV) 

1 In (T/T K ) 

hi {T /Tk) 
\ 1/2 

Tk 



TABLE II. Low temperature behavior of the impurity 'orbital susceptibility' Xim P = 
d Fimp/d/S. , specific heat Ci mp and resistivity i?i mp of a TLS. The energy scales, T* , 
T, A, are assumed to be smaller than the Kondo energy Tk. 
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VI. GENERALIZATION TO THE CASE OF AN M- STATE SYSTEM 



In the previous Section we discussed the special case of two-level systems. We 
have seen that at small energy scales these can be discribed by the two-channel 
Kondo model. The, TLS model is appropriate to describe tunneling centers in 
amorphous metals,E3 where the probability of having a three-state system is much 
smaller than that of having a TLS. However, it breaks down in systems, where 
the tunneling centers are formed by some substitutional impurities, and the heavy 
particle is tunneling between 3, 6 or 8 equivalent positions. Tipical examples of 
such aikttsa are the narrow gap semiconductor Pbi-^Ge^Te or insulating Ki^Li^d 
alloys e3 c3 Three-state systems have also been observed in MOSFET devices where 
most of the tunneling defects are TLS's due to some random structured Therefore 
it is natural to ask what is the low temperature behavior of an M-level system 
(MLS) which strongly interacts with the conduction electrons. 

To answer this question we first have to generalize our TLS model. Assuming that 
the temperature (or the relevant energy scale) is low enough and thus the motion 
of the heavy particle is now restricted to the lowest lying M states co rres ponding 



to the M spatial positions of the heavy particle the TLS Hamiltonian (2.1) will be 
replaced by 

M 

H hp = ]T b+AVbj , (6.1) 

where bf creates a heavy particle at site i and A 4 - 7 is the tunneling amplitude 
between positions i and j. If we assume that no external stress is present and 
that the M positions are completely equivalent thus A M = 0. This assumption is, 
however not necessary for the considerations below. 

Similarlly to the TLS case the we consider a general electron spin s taking the 
values s = 1 . . . Nf and we shall develop a 1/Nf expansion around the fixed point. 
The kinetic part of the Hamiltonian remains unchanged and the most general two- 
particle interaction generated by some effective interaction (screened Coulomb in- 
teraction or a pseudopotential) between the heavy particle and the conduction 
electrons now takes the form: 

flint = b i V nil b j a tns a e'm S ■ (6.2) 

€,€ 

For the sake of simplicity we also assume a constant density of states per flavor 
go between the high- and low-energy cutoffs D and —D, for all flavor numbers. 
Naturally, both the couplings V^ m and the tunneling amplitudes A 1 -? are connected 
by the symmetry properties of the MLS which will be exploited later on. 

In the MLS case the pseudofermion propagator Q l i(uS) and the vertex function 
r lJ (cj) are M x M matrices in the pseudofermion indices. Calculating these func- 
tions in a perturbative way it turns out that th ey d o not satisfy the simple multi- 
plicative renormalization group (RG) equations (|4~l| — 14.3|) and the following gener- 
alized renormalization group transformation must be used (T = 0) £U 

G{uj, v', A', D') = A G{u, v, A, D) A+ , 

T(u,v', A',D') = [A + }~ 1 T(uj,v, A, D) A' 1 , (6.3) 

where the matrix notations QoV^ n — > v 1 ^ — > v, ^mn ~ * E lJ ~~ ¥ —> an< ^ ^ ~ * ^ have 
been introduced, D' stands for the scaled bandwidth and A is an M x M matrix 
acting in the heavy particle indices. Note that A = A', D' /D) is independent 
of the dynamical variable ui. 

This transformation has the following properties: 

(i) For A % 3 = Z 1 > 2 S l; > it is identical with the usual multiplicative renormalization 
group transformation. 



(ii) The symmetry Eq. (6.3) is generated by the 'microscopical' transformation 
v -> {A+^vA- 1 , A~U (A+)- 1 A J 4- 1 , g Q -> Ag A+, g a being the bare 
heavy particle Green's function. 
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(iii) The transformations (J6 .3|) form a semi-group 



(iv) Eq. leaves the Hamiltonian Hermitian. 



While for finite D/D' the matrix A has a rather complicated structure, still there 
exist such physical quantities like the electronic scatteri ng r ate or the impurity free 



energy, which are invariant under the transformation (3.3) and can be calculated 



easily. It is important to note that for an infinitesimal change of D the matrix A 



can be chosen to be Hermitian and Eq. (6.3) can be cast in the form of a scaling 
equation for the dimcnsionlcss couplings v 10 . 

In the following we assume for the sake of simplicity that the relevant energy 
variable is u, i.e., lu 3> |A y |,T. In this case the inverse heavy particle propagator 
and the vertex functions can be expressed in the next to leading logarithmic ordenla 
as 

(g- 1 )^ = uj ft - A ij + N f In- (S ij u Tr{v kl v lk } - Tr{v ik A kl v lj }) 

go ri = v ij - In — ([v ik ,v kj ] - N f Tr{v ik y! j }v kl ^ , (6.4) 

where [ , ] denotes the commutator. Plugging ( |6.4| ) into Eq. ( |6.3[ ) one can easily 
generate the scaling equations in a selfconsistent way and after some algebra one 
finds: 



dA» 1 

dv ij 
dx 



\v ik ,v kj ] + 



Tr{v ik v kl }A lj + A ik Tr{v kl v lj } - 2Tr{v ik A kl v! j } 

^N f (2Ti{v ik v l3 }v kl - Ti{v ik v kl }v lj - v ik Tr{v kl v lj }^ , (6.5) 



where x = In (Dq/D) denotes the scaling variable. It is important to note that in the 
next to leading logarithmic level the scaling of the 's is completely independent 
of the splittings A y , while the scaling of these latters is driven by the couplings 

These scaling equations are very complicated and, apart from some particular 
cases, one can solve them only numerically. However, to exploit the symmetry 
properties of the MLS it is useful first to introduce a site representation in the orbital 
indices of the conduction electrons. We do this by taking some linear combinations 
of the most strongly scattered angular momentum channels and hybridize them 
using group theoretical methods. For a regular 3-state system in the xy plane and 
a free electron band, e.g., one can use the three orthogonal electron states: 

|e2> = 73 |es)_ 76 |ePx) + 72 |ePy> ' (6 ' 6) 

\e3) = ±\es)--L\e Px )--L\ep y ), 

where the states \p x ) and \p y ) are defined in the usual way from the angular mo- 
mentum states (r|e I m) = /(d 3 k/(27r) 3 )y im (k)e lkr <5(e - e k ) with I = 1. Working 
with only those electron states which are directed to the impurity positions the 
v becomes an M 4 -dimensional tensor. However, the number of independent cou- 
plings is largely reduced by symmetry. For a heavy particle tunneling between the 
six corners of a regular octahedron, e.g., the 6 4 = 1296 couplings may be replaced 
by 32 independent couplings which makes a numerical solution reasonably fast. 
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FIG. 12. Scaling of the norm of the dimensionless couplings, u = (dashed 
line), and of the algebra coefficient a (continuous line) for a 6-state system with Nf = 2. 
The initial couplings have been chosen to be v{{ = 0.8, W22 — 0-2, «66 = 0-1) v i2 = 0.05, 
v\l = 0.03, W21 = v\l = 0.0005, v\l = 0.005, v\t = 0.003, v\l = v\l = 0.0005. The other 
nonzero couplings have been generated by symmetry transformations. 

In Fig. Owe show the typical scaling of the norm of the dimensionless couplings, 
u = S ll w (dashed line), where the initial couplings have been estimated simi- 



larly to Refs. |18|, |20| One can show similarly to the multichannel Kondo problem 
that both the infinite and the weak coupling fixed points .ape. unstable and the sys- 
tem scales to an intermediate strong coupling fixed point .EjEj'lj The Kondo energy 
can be defined as the crossover energy from the weak to strong coupling regimes: 
Tk = D e~ Xc , x c being the crossover value of the scaling parameter. We foud that 
for realistic initial parameters this Kondo temperature can easily be found in the 
experimentally observable region, Tk ~ 1 — 10 if. 

Similarly to the TLS case we would like to determine the properties of the MLS 
model in the regime T* < T, u <C Tk- (Below therfceezing temperature T* the 
dynamics of the MLS becomes trivial in most cases.E3) To that purpose we have 
to determine the low-temperature effective model, i.e., the stable fixed points of 



Eq. (5.1). Obviously, the scaling equation has an infinite number of fixed points. 
Similarly to the TLS case any representation of the SU (M) Lie algebra will define 
a possible fixed point. Trivially, the TLS fixed points, where two MLS states are 



completely decoupled from the other M — 2 states, are fixed points of Eq. ( |6.5D as 
well. Numerical investigations show that, apart from the fixed point where the v^'s 
realize the defining representation of the SU{M) Lie algebra, all these fixed points 
are unstable. The special fixed point corresponding to the definingjjepresentation 
is equivalent to the the SU{M) x SU(Nf) Coqblin-Scrieffer modelO 

In the following we shall show that the above-mentioned Coqblin-Schrieffer model 
is really stable. To this end we first remark that the operators 6^ J2k Unm are 
invariant under scaling. Therefore the u y 's can be written as = + M^ n , 
where the matrix M is a constant of motion depending on the initial parameters 



and Yli^nm — 0- Then one can easily show that the right-hand side of Eq. (6.5) 
disappears provided 



J_ / iji 
N f \ 



(6.7) 



where the L y 's are unitary equivalent to the generators of the SU(M) Lie algebra, 

[L ij ,L M ] =8 a L kl -5 k3 L a . (6.8) 

with L^ m ~ 5 l m 5 3 n — jj8 % i 5 nm . To demonstrate that an MLS scales to this fixed point 
in Fig. |l2| we show the scaling of the 'algebra coefficient' a — J^i j k i\ > v kl ] — 
NfS v i + N fS iy \\, which measures in a natural way how well the fixed point 



algebra (6.8) is satisfied. For D < Tk (i.e. for x > x c ) the coefficient a scales to 
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zero. Thus we conclude that below the Kondo temperature an MLS scales to the 
SU(M) x SU(Nf) Coqblin-Schrieffer model, which is a non-Fermi4imtid model and 
has a different scaling behavio r th an the 2-channel Kondo model.cJE3 

To show that fixed point (|6.7|) is stable and to analyze its operator content 
we follow the lines of the stability analysis of the TLS fixed point and write the 
deviations of the couplings from their fixed point value in form 

*L U = ( JL £ ) , (6-9) 



(t ji )+ K 13 

where the couplings g* 3 , f y , and g %3 are M X M, M X oo, and oo x oo matrices, 
respectively. Like the TLS case the linearized equations for g 13 , t 13 , and g %3 decouple 
completely, 

JL = _ _ M /) , (6.10) 



dg il 
dx 



= — — ([L ik , g kl ] + [g ik ,L kl ]) + — \2S u g kk + 2L jk Tr{g i3 L kl + L ij g kl } 
- 2Mg il - L i3 Tr{g> k L kl + L jk g kl } - Tr{g i3 L jk + L ij g 3k }L kl \ (6.11) 

rjfH 1 1 

— = (L ik t kl - L w t tt ) + — (5 il t kk - Mt a ) , (6.12) 

and they can be solved exactly due to the simple structure of the L^'s. The 
linearized equations have an infinite number of zero modes; a finite number of 
them correspond to potential scattering while the others can be identified with 
the generators of the unitary transformations connecting the different possible M- 
dimcnsional subspaces where the SU(M) Lie-algebra is realized. These 0-modes 



can be shown, of course, to leave the Lie-algebra (6.8) invariant. All the other 



modes can be shown to be irrelevant, which proves the stability of the fixed point 



(6.7) 



As we have pointed out in the previous Section, the universal properties of a 
model are determined by the operator content of its stable fixed points. Therefore, 
in order to determine the properties of the MLS in the region T* < T < Tk, 
T* being the freezing temp erat ure, we have to deter mine the operator content 
associated to the fixed point (6.7). The analysis of Eq. J6.12D shows that for M > 3 



the leading irrelevant operators can be written as 

Oi - ( {c l )+ °q ) , (6.13) 

where the C ij 's satisfy (C k L - C%£) = 0. These operators scale like - D Al with 
A; = (M — 1)/Nf, and they describe scattering between channels which are not 
taken into account in the usual Coqblin-Schrieffer model. While they dominate the 
thermodynamical quantities like the specific heat, e.g., which scales as Ci mp ~ T 2Al 
they do not contribute to the resistivity, which scales like ~ T Ab1 with A s ; = M/Nf, 
and is determined by subleading operators of the form 

Osi - ( Q q ^ ) , (6.14) 

where the matrices Q lj and satisf y Y", Q u — S m — and = Q% m . It is 
important to note that the operators ( |6 . 13 ) do not exist in the TLS case (M = 2), 



which explains why the last-energy properties of a TLS can be described by the 
two-channel Kondo model.E3 

It is an open qestion whether the leading irrelevant operators found in the next to 
leading logarithmic approximation really exsist for the physical case Nf = 2. Our 
results become, like the TL^case, exact in the Nf — ► oo limit. However, wejfixpect 
from conformal field theoryEil and noncrossing approximation (NCA) resultsEa that 
the 1/Nf expansion breaks down at Nf = M . Comparing the previous results with 
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the exponents of the conformal field theory and recent Bethe Ansatz calculations^ 
SU(M) x SU(Nj) model we can identify the exponents M/N } and (M - 1)/N f 
with the exact ones M/(M + N f ) and (M - 1)/(M + N f ). The expansion of the 
latters with respect to 1/Nf breaks down at M — Nf. Similarly, the NCA gives 
the result that the M > Nf and M < Nf models behave in a different way and 
that the M = Nf case is marginal. We expect that, similarly to the TLS case 
or the anisotropic Kondo model, the fixed point symmetry remains the same even 
for Nf < M , however, the dimension of the different operators and the operator 
content might drastically change as we go to the region M > Nf. Therefore the 
role and the physical interpretation of the new leading irrelevant operators in the 
Nf = 2 case is still rather obscur and needs some clarification. 

Up to now we have ignored the effect of the MLS splitting. These splittings, as 
mentioned before, can be taken into account as infraread cutoffs and they finally 
stop the scaling towards the SU(M) x SU(Nf) non-Fermi liquid fixed point at 
the freezing temperature now defined by T* ~ max i t3 - \A V (T*)\. Much below 
this scale the MLS is frozen into its ground state. Usually this ground state is 
nondegenerateL3 and therefor e, if T* > Tk then the non-Fermi liquid properties 
associated to the fixed point (6.7) can not be observed. Therefore, the estimation 
of the splitting parameters and their renormalization under the scaling procedure 
is very important. 
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FIG. 13. Scaling of the dimensionless hopping amplitude, A 12 /Do for the same 6-state 
system as in Fig. 1. Inset: Numbering of the sites of the 6-state system. 

To estimate the renormalization of A y we solved Eq. ( |6.5D numerically. As 
shown in Fig. [l^, for realistic model parameters A 12 (Tk )/A (Do) can be as small 
as ~ 1CP 3 , and therefore, even for very large splittings A 12 ~ 10QK the splitting 
is strongly reduced and since A 12 (Tk) <C Tk the MLS can get easily into the 
vicinity of the 2-channel Coqblin-Schricffcr fixed point. We stress at this point that, 
similarly to the TLS asymmetry energy A z , the renormalization of the diagonal 
part of the pseudofermion (MLS) self-energy ~ (A™ — A-" ) is much smaller than 
that of the tunneling amplitude A y (i ^ j). However, for an MLS, (if there is 
no stress in the material, and the concentration of MLS's is small), the previous 
ones may be quite small due to the symmetry of the MLS, and the formation of 
a Kondo- like ground state with T* < Tk, i.e., with observable non- Fermi liquid 
behavior becomes quite reasonable. For a TLS in an amourphous material no such 
symmetry exsists and the non- Fermi liquid properties are very probably smeard out 
in the bulk system unless a very peculiar distribution of the TLS splittings exists. 

Unfortunately, metallic samples containing such highly symmetrical MLS's due 
to some interstitials are not too numerous. One example of them displaying a more 
or less unambiffous Kondo effect under preassure is the narrow gap semiconductor 
Pbi^^Ge^TeH In this material the relatively small Ge 2+ ions form 8- level systems 
and interact with the conduction electrons. Alas, this material is a very complicated 
one and its preparation was, at the time of the measurements, not well controlled 
neither. In these samples the concentration of the MLS's is rather high and their 
interaction probably results in the nonvanishing of the difference of the MLS site 
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energies, A M — A 3 - 7 , which hinders the observation of the non- Fermi liquid properties. 
Moreover in the Pbi_ x Ge x Te alloy the spin-orbit interaction is also large and may 
result in a cross-scattering between the two spin- channels. This process also spoils 
the non-Fermi liquid behavior. Thus, in order to observe the predicted non-Fermi 
liquid property of an MLS one should find some better candidate. 

One of the other possibilities is to use nanotechnology to construct devices to 
realize the SU(M) x SU{2) models. Such a realization of the SU(3) x SU(2) model 
has been suggested in Ref. [53] by means of a double dot structure. 

VII. SHORT REVIEW OF NEW THEORETICAL DEVELOPMENTS 

In the previous Sections we mainly concentrated on the mapping of an MLS 
model to the SU(M) x SU {Nf) Coqblin-Schrieffer model by using the multiplicative 
renormalization group technique. These issues, however, constitute only a fraction 
of the presently available results for tunneling Kondo impurities. The aim of the 
present Section is to review in a concise way such techniques and results that lie 
out of the scope of the present work. To get a deeper understanding of them the 
interested reader is referred to the original works. 

A. Path integral treatment 

One of the drawbacks of the multiplicative renormalization group treatment of 
the TLS (MLS) model is that it is only reliable in the small coupling region where 
the oerturbative expansion is meaningful. On the other hand, simple estimations 
showta that for reasonable TLS parameters the screening interaction v z might be 
quite strong, v z ~ 0.4, where the multiplicative renormalization group is already 
not applicable. This discrepancy of the multiplicative renormalization group, has 
partially been cured in Ref. |2g where a Yuval- Anderson type scaling analysisEj has 
been applied to the TLS model. 

In this approach, as a first step, one devides the Hamiltonian into two parts, 
-ffnoflip and -ffflip, where all the terms associated to the flip of the TLS (v x , v v , 
A x , and A v ) are grouped in the part -ffflip. The main observation is that -ffnofiip 
contains all the large parameters of the model and can be solved exactly. Therefore, 
calculating the partition function Z of the model, one can develop a perturbation 
theory in the small TLS-flipping interactions. Due to the relatively simple structure 
of iJnoflip one can integrate out the electronic degrees of freedom in each term of 
the perturbation series. As a result of this procedure one is left with an effective 
action for the TLS pseudospin only. 

As usually, this effective action is fairly complicated and containes logarithmic 
retarded interactions between the different TLS flips. However, one can use it to 
generate renormalization group equations even if v z is large, by requiring the in- 
variance of the partition function under the renormalization group transformation. 
The prcedure is, however, very sensitive to the exact/way of elimination of the high 
energy degrees of freedom and needs lots of care.EJ Several different elimination 
shemes exist and some of them leave the partition function invariant but not the 
different dynamical Green's functions. Another disadvantage is that the expansion 
of Z is not systematic and by construction does not treat the different couplings on 
equal footing. Up to now no-one succeded to derive the next to leading logarithmic 
scaling equations with this method. 

The TLS model has also been studied recently by a bosonization technique com, 
bined with scaling analysis with somewhat different results from the previous ones£3 
In this method the electronic degrees of freedom are represented by their bosonic 
charge (spin, flavor, and spin-flavor) density fluctuations which makes the problem 
easier to handle. Very probably most of the aformentioned differences arize from 
the fact that in Refs. [58] such an elimination scheme has been used that mixes to, 
gether the TLS and conduction electron degrees of freedom in an artificial wayE3 
On the other hand, the bosonization procedure uses a very special cutoff scheme 
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and it is presently not clear for us if the anomalous dimension of the different cou- 
plings is correctly given by it at a finite phase shift for an arbitrary process such as 
multi-electron scattering. 

The Anderson- Yuval technique has also been SHCcesfully applied to the problem 
of the tunneling of a heavy particle on a lattice.Ea For the noncommutative case, 
however, the scaling equations become very involved and little physical informa- 
tion can be obtained from them. This problem has also been studied using the 
multiplicative renormalization group method in Ref. |oC| . 



B. Role of the excited states 



In the previous Sections we simply ignored the excited states of the TLS by the 
handwaving argument that they are completely frozen out at the low temperatures 
we are interested in and only the lowest-lying nearly degenerate states must be kept. 
These excited states, however, influence the low-energy behavior of the TLS. Even if 
the temperature is low virtual processes consisting of several 'assisted hoppings' to 
the excited states like the one shown in Fig. |l4| are allowed, generating an effective 
assisted tunneling process from one well to the other. Since the overlap of the states 
|Z) (|r)) and \ex) is large compared to that of \l) and |r), these processes cannot be 
neglected. To take them into account one can easily generalize the Hamiltonians 



( pTlD and (gj) to 

#TLS = J2 Eib t b i > 

i 

H int = J2 ^X'b+bja^u , (7-1) 



e,n,e ,n 



where now the operators bf create a state with the heavy particle in the exact 
eigenstate of the double well potential \i) with energy J5{, where E\ = s» E 2 <C 
£/3 < E4, < ... 



U ( r ) 




FIG. 14. Fig. a: Sketch of a double well potential and the states \r >, \l >, and the 
first excited state \ex >. Fig. b: A diagram generating the assisted tunneling via virtual 
transitions to the excited states. 



A careful estimation of the assisted hopping matrix elements and a multistep 
scaling analysisEl reveal that for a typical TLS the presence of the excited states may 
increase the Kondo temperature by a factor of 10—100 thereby giving Tk ~ 1 — 10K 
which is in much better agreement with the experimental results as our previous 
estimation (4.14). It has also been shown in Ref. |2(] that at low temperatures the 
assisted hopping to the excit ed st ates can be taken into account by the redefinition 
of the couplings VJL, in Eq. (2.4). Thus one can ignore formally the excited states 
but then the TLS-electron couplings vf 1 should be viewed as some effective couplings 
incorporating the influance of them as well. 
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C. A TLS with spin 



To obtain a non-Fcrmi-liquid scaling of the TLS it has been crucial that the 
scattering of the conduction electrons be diagonal in ther real spin (~ flavor), s — 
1, ..,Nf. In this case the Hamiltonian possessed an additional SU(Nf) symmetry 
associated to the flavor degrees of freedom. From the group theoretical point of 
viewcfl it is this additional SU(Nf) symmetry that is responsible for the non- Fermi 
liquid behavior of the model. Therefore it is natural to ask what happens if the 
tunneling particle has also a spin which is coupled to the conduction electrons by 
a local exchange interaction ~ JS l a l (R), cr l (R) being the local spin density of the 
conduction electrons at th e heavy particle's location R. In this case, in addition to 
the usual interaction (|2.4|) an exchange interaction term appears of the form 



E 



J nn' u n 



+ % 



>„>). 



(7.2) 



i,f-i,s,s ,n,n 



This interaction obviously breaks the SU (2) symmetry in the electronic spin, and 
it mixes together the two electronic spin channels. A simpler version of Eq. (7.2) 
has already been studid in Ref. 25 in the context of magnetic impurities in crystal 

showalU 



(7.2 



A straightforward renormalization group study of the model Eq. 
that, as one could guess it from the outset, the coupling to the impurity spin always 
drives the system to a Fermi liquid ground state even in the absence of the TLS 
splitting A. The nature of this ground state may, however, depend on the initial 
parameters of the model. 
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Fig. Ida 



Fig. M 



FIG. 15. The scaling trajectory of the relative norms ||u||/(||j|| + IMI) (dashed line) 
and ||i||/(||j|| + IMI) (solid line) in a logarithmic scale for two different set of initial 
couplings. The scaling trajectories go from right to left, and x c = \xi{Do/Tk)- In Fig. a 
the exchange interaction dominates the fixed point behavior of the model. For the small 
exchange coupling case in Fig.b. a different fixed point appears, characterized by a finite 
universal ratio ||i||/|M|- 



In Fig. |l5j the scaling of the relative norms of the dimensionless couplings 
and vj 1 is shown for two sets of initial couplings in the leading logarithmic approx- 
imation. In Fig. |l5] the model scales to a fixed point where the couplings can 
be neglected compared to the exchange couplings. In this case the generalized ex- 
change coupling ([7j) generates a dynamical splitting for the TLS and therefore the 
orbital motion of the TLS is frozen out at low temperatures, and the low-energy 
behavior of the model can be appropriately described by a spin 1/2 single channel 
magnetic Kondo model.til 

On the other hand in Fig. |l5| b we show another situation where the ratio ||i||/||f || 
scales to a finite value. In this case at the fixed point both the orbital and magnetic 
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degrees of freedom are coupled to the TLS in a symmetrical way and the low-energy 
effective model can be shown to be equivalent to the single channel SU (4) Coqblin- 
Schrieffer model. 

The tunneling spin problem has also been treated very recently usinsj-conformal 
field theory methods which resulted similar results to the present ones.Ea 



VIII. REVIEW OF THE EXPERIMENTAL SITUATION 



Presently there are only a few systems where Kondo-like anomalies attributed to 
TLS's or other tunneling centers have been observed. 

Logarithmic anomalies due to TLS's formed probably by dislocations have been 
found in some amorphous systems like Ali-^Mg^ or Ali-^GeJ^j In these materials 
the typical splitting of the TLS's is large (~ 1 — 10K) and a wide distribution of 
the TLS splitting is present, therefore the non-Fermi-liquid regime can never be 
observed. These anomalies are not sensitive to magnetic field and they disappear 
upon annealing. 

Another extensively studied system containing Kondo-like tunneling centers is 
the alloy Pbi_ x Ge x Te that we already discussed in detail at the end of Section [V^ . 
While in some circumstances_an unambigous logarithmic anomaly due to TLS's can 
be observed in this materialJEj it is too complex from the theoretical point of view 
and the interplay of different fenomena (phonon softening, resonant scattering from 
phonons, ferroelectric phase transitions) make it very difficult to analyse. 



FIG. 16. Universal scaling behavior of the unannealed Cu point contact spectrum from 
Ref. |I| 

The most promising experiments co ncern ing TLS Kondo impurities have been 
performed on metallic point contacts qc3 In these experiments only scattering 
on a few TLS's in the contact region influences the signal of the device, and the 
differential conductance of the point contact is directly proportional to the tem-i 
perature and energy dependent scattering rate of the electrons on these TLS's,E-i 
G(V, T) = J^r ~ t(tv) • ^ n ^ese measurements not only a logarithmic temperature 
and voltage dependence at high energies has been observed but also a universal scal- 
ing ~ [max(T, V^)] 1 / 2 of the conductance curves shown in Fig. ^| has been found, 
The two-channel Kondo scaling curves calculated by noncrossing approximation^ 
and conformal field theorjCa fit very well the observed scaling curves. 

Altshuler et al. suggested that a similar scalieg behavior could result from 
electron-electron interaction in localization theory.t3 However this explanation has 
several constituents which are clearly in contradiction with the experimental facts 
(an anomalously short estimated mean free path being in contradiction with the 
good resolution of the phonon peak, disappearance of the signal with increased 
static disorder, magnetic field independence of the signal in some experimentsES) Ea 
We should mention, however, that while one could see in some experiments the 
cutoff of the scaling by soni|e-|intrinsic splitting at low energies in complete agree- 
ment with the TLS theory,E3 in others no violation of the universal scaling has 
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been observed at T ~ IK. This should imply that TLS's with small renormalized 
splittings A* < IK are very likely to occur in the latter systems. This might be 
due to some special origin of the TLS's but it is presently not well-understood. 




V [mV] 



FIG. 17. The squares and triangles correspond to AG(V) calculated from curve 3 in 
Fig. 2 and the curve of Fig. 4 of Ref. j7^. The continuous curves are fits with the Kondo the- 
ory (the dashed line has been fitted by scattering from two TLS's). The inset shows curves 
calculated with the Kozub-Kulik theoryO incompatible with the experimental results. 

Finally, we discuss the results of a very recent measurement on amorphous break 
junctions by Keijsers et al£2l giving a very nice possibility to observe directly a 
single noncommutativc TLS. In these experiments Keijsers et al. found that the 
zero-bias anomaly G(V) (R(V)) of the point contact fluctuated slowly between two 
values, G\(V) and G^[V). Since the zero-bias anomaly is thought to be due to 
the fast two-level systems present in the contact area these fluctuations can be 
explained in terms of a slow fluctuator that influances the parameters of the fast 
TLS'sJia One can also argue that the difference AG = Gi{V) - G 2 {VX should be 
due to only a few possibly only one fast TLS with modified parameters O Since it is 
the asymmetry parameter A 2 of a TLS that is most sensitive to external stress, it 
should be possible to fit the difference AGlK) by changing the splitting parameters 
A; -> A^ for just a few TLS Kondo curvesO Our fit to AG(V) is shown in Fig. [l7|, 
and is in very good agreement with the experimental curves, thus supporting the 
existence of fast Kondo-like two-level systems in amorphous metals. 

IX. CONCLUSIONS 

In the present paper we tried to give a comprehensive and somewhat detailed 
review of the present understanding of TLS Kondo impurities. While we mainly 
concentrated on the multiplicative renormalization group approach and to the map- 
ping of M-level systems to the much simpler SU(M) x SU(Nf) model we also at- 
tempted to review shortly other theoretical and experimental results in the field. 
To understand these more in detail the interested reader is adviced to consult the 
given References directly. Of course, this review is far from being complete since 
the TLS literature is growing very fast. The role of excited states is still not well 
understood, very little is known about the influence of Kondo TLS's on the other 
properties of the material and the properties of M-level systems. 

We hope that we managed to convince the reader that tunneling Kondo impurities 
deserve further studies and provide very promising realizations for the two-channel 
Kondo model, even though the discussions are not settled about them yet. 
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